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■ Abstract. We follow the evolution of the Black Hole mass function in 
\ the assumption that mass accretion onto a massive BH is the powering 

J ^ ■ mechanism of AGN activity. With reasonable assumptions suggested by 

O , current knowledge, it is possible to reproduce the Black Hole mass func- 

' tion derived from local bulges. In particular it is found that the majority 

I , of the mass in relic Black Holes was produced during AGN activity and 

ff^ ' that the ratio between obscured and unobscured quasars is between 1 and 

2, much lower than required by synthesis models of the X-ray background. 
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^ ' 1. Introduction 

o : 

^~~] ' It has long been suspected that the most luminous Active Galactic Nuclei (AGN) 

, are powered by accretion of matter onto massive black holes (BH). This belief, 

^ ' combined with the observed evolution of the space-density of AGN suggests that 

p . massive BHs should be present in local galaxies as relics of past activity. Indeed, 

I ' the current observational evidence indicates that most, possibly all, luminous 

2 . galaxies host a massive BH at their centres and its mass correlates with the 

\ luminosity and stellar velocity dispersion of the bulge in which it resides (e.g. 

! Richstone et al. 1998, Ferrarese & Merritt 2000). Soltan (1982) and Chokshi 

\ & Turner (1992) derived the local mass density in relic BHs, /jbh > expected 

• j . from past activity using optical number counts of AGNs. More recently, Fabian 

^ \ Sz Iwasawa (1999) estimated pbh from the energy density of the hard X-Ray 

■ Background (XRB) and this value is in agreement with the observed local BH 
density derived by Merritt &; Ferrarese (2001) using the BH-mass vs Bulge-mass 
relation. However it seems that the AGNs found in optical and soft X-ray surveys 
are not enough to account for the local pBU and this fact prompted Fabian 
& Iwasawa to talk about "obscured accretion", i.e. obscured AGN activity. 
The existence of obscured AGNs (the so-called "type 2" objects) is required 
to account for the XRB spectrum with integrated emission from AGNs (e.g. 
Gilli, Salvati & Hasinger 2001 and references therein): the ratio between type 
2 and unobscured type 1 objects (hereafter R2-1) is usually in the range 4-10, 
implying that both optical and Xray surveys are missing a large fraction of 
the AGN population. With a completely independent argument, the missing 
Lyman edge in blue quasars implies that R2-1 is at least > 0.5 (Maiolino et 
al. 2001). Locally, low luminosity obscured AGNs are well known (Seyfert 2 
galaxies) and in the ratio of about 4:1 with respect to type 1 objects (Maiolino 
& Rieke 1995). However, at higher redshifts few type 2 quasars (high luminosity 
objects) are known, a number which is largely insufficient to account for the 
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Figure 1. (a) AGN luminosity function at z = 1.0 (solid line) and 
z = 2.5 (dashed line). This is obtained by combining the XLF by Gilli, 
Salvati & Hasinger 2001 (low L) and the OLF by Boyle et al. 2000 (high 
L). (b) Local BH Mass Function (solid line). The uncertainties are 
represented by the gray area. The dashed line is the BHMF obtained 
assuming that ax = 0- 



spectral shape of the XRB. What is the relationship between type 2 AGNs 
and massive black holes? The arguments based on pbh indicates that some 
obscured AGNs must exist in order to account for the fossil mass but nothing can 
be told about the luminosity of these obscured objects: are many low luminosity 
AGNs needed or could fewer high luminosity quasars suffice? Is it really true 
that AGN activity can explain the local BHs, not only globally, but also in 
terms of their mass distribution? This paper will try to answer these questions 
by comparing the local BH mass function (BHMF) expected from AGN activity 
with the local BHMF inferred from galaxy bulges. Throughout this paper we 
assume = 70km/s/Mpc~^, f^u = 1 and JIa = 0. 



2. The continuity equation 

We consider the continuity equation which regulates the time evolution of the 
BH mass function with the Small & Blandford (1992) formalism. If N{M,t) 
is the BH mass function (N{M, t)dM is the number of BHs per unit comoving 
volume (Mpc"'^) at cosmic time t), the continuity equation can be written as: 



dN{M, t) d 



dt dM 



(^N{M, t) < M{M, t) >) = S{M, t) (1) 



where < M{M,t) > is the "average" accretion rate and S{M,t) is the source 
function. As in Small &; Blandford (1992) we assume that a massive BH is 
"born" when its mass is M > Mmin and that no BH is created or destroyed (i.e. 
we neglect any merging process) . With these assumptions the source function is 
simply S{M, t) = So6{M - Mmin) where 6 is the Dirac delta function. If AGNs 
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are powered by accretion onto BHs, we can relate the AGN Luminosity Function 
4){L,t) ((/>(L, t)dlnL is the number of AGNs per unit comoving volume (Mpc~^) 
at cosmic time t) to the BH mass function 

t)d \nL = 5{M, t)N{M, t)dM (2) 

where 5{M,t) is the fraction of BHs with mass M active at time t. If a BH is 
accreting at a fraction A of the Eddington rate its emitted luminosity is then 
L = XMc^/tE = eMc^ where tE is the Eddington time, e is the accretion 
efficiency and M is the accretion rate. Since < M{M,t) >= 6{M,t)M{M,t), 
we can write: 

N{M, t) < M{M, t) >= [^{L, t)]^^^M^ (3) 



dN{M,t) _ c^A^ 



dt £t\ 



£tE ^ " tE 

d(t>{L,t) 



dL 



(4) 



E 



which can be easily integrated given the AGN luminosity function AND the 
initial conditions. When A is dependent on M the equation is more complex but 
still easily solvable. As initial conditions, we assume that at a given redshift Zo 
(to = t{zo)) ALL Black Holes are active, i.e. that 5{M,t{zo)) = 1. The only 
constraint is that 6{M,t) = (j){L,t)/MN{M,t) < 1. 

A fundamental ingredient is the AGN luminosity function (LF) expressed 
in terms of the total AGN luminosity, L. We first consider the LF of optically 
selected quasars (OLF in the B band; Boyle et al. 2000) and the LF of soft X- 
ray selected AGNs (XLF in 0.5-2 keV band; Gilli et al. 2001, it represents only 
unabsorbed AGNs and is derived using the XLF by Miyaji, Hasinger & Schmidt 
2000). Both LFs consider only type 1 (unobscured) AGNs. In principle the LF 
in one spectral band should be enough to obtain 4>{L, t) but in this paper we 
consider more instructive to combine the above LFs using Boyle's LF to describe 
high L AGNs and Gilli's one to describe low L AGNs. A more refined treatment 
will be presented elsewhere. The AGN luminosity function for type 1 objects is 
plotted in Fig. ||a at two redshift values. We have assumed L/Lx = 30, for low 
L objects and L/vqLi/-q = 15 as typical of quasars (e.g. Elvis et al. 1994). 



3. The local BH mass function 

To estimate the local BH mass function we follow the method by Salucci et al. 
(1999). The BHMF is derived from the Galaxy LF through the following steps: 
(1) we consider the LF of galaxies in the B band divided by morphological type 
(Folkes et al. 1999); (2) for each morphological type we derive the LF of bulges 
applying the bulge/disk ratio by Simien & de Vaucouleurs (1985); (3) we obtain 
the mass function of bulges using the mass-to-light ratio by Fukugita, Hogan, 
& Peebles (1998); (4) we convolve the bulge mass function with the probability 
distribution of the ratio MBH/-^Buige to obtain the BH mass function. An 
important assumption in the last point is that the observed scatter in the Mbh- 
-^Buige relation is real and not due to observational uncertainties. We then 
assume a gaussian probability distribution P(logx) (x = MBH/-^Buigc) with 
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Figure 2. (a) Local BH mass function from galaxy bulges (solid line 
in grey area) compared with BH mass function obtained considering 
only type 1 AGNs (solid line). The dotted lines are the contributions 
from low and high L AGNs. The dot-dashed line is the BH mass 
function at the redshift Zo=3 when all BHs were active. The dashed 
line represents the BHMF from AGN activity assuming that only 10% 
of the BHs were active at Zo=5. (b)The contribution from high L AGNs 
is computed for R2_i,b=0 (solid line), 1 (dotted line), 2 (dashed line). 




Figure 3. (a) BHMF from AGN activity computed with R2_i,x=4 
for L/Lx = 30 (solid line), 20 (short-dashed line), 10 (long-dashed line). 
The dotted line represent the model with R2_i,x=0. (b) The dotted 
and solid lines represent models with L/Lx = 30, L/u-qLu^ = 15 and 
R2-i,x=4. In the model represented by the dotted line all AGNs are 
emitting at their Eddington luminosity, while with the solid line AGNs 
are emitting at a fraction of their Eddington luminosity which depends 
on M as specified in the figure. For this model R2-i,b=1- 
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< logx >= —2.93 and ax = 0.45 as derived by Merritt & Ferrarese (2001). The 
BHMF is presented in Fig. ^3 and the "grey" region represent an error estimate 
on the BH mass function obtained by varying the parameters used in each of the 
above steps. The BHMF has pbh = (5 ± 2) x W^Mq Mpc~^ and this value is in 
agreement with pBH = (3 -i- 6) x IO^Mq Mpc~^ obtained by Fabian & Iwasawa 
(1999) and Salucci et al. (1999) from the hard XRB. 

4. Model Results 

Model results are shown in Fig. ^a. The solid line in the grey area is the BH 
mass function derived from the galaxy mass function as indicated in the previous 
section. The solid line is the BH mass function expected from type 1 AGNs only 
with a starting redshift of Zo = 3. The starting redshift at which we make the 
assumption of 5{M,to) = 1 (or alternatively the redshift at which we select the 
active BHs whose evolution we will follow) is not important for the final results. 
Indeed the dotted line representing the BHMF obtained assuming Zo=5 and 
6 = 0.1, is not significantly different from the one with Zo=3 and 5 = 1. The 
only constraint is that Zq > z^ax, the redshift of maximum AGN activity. From 
Fig. |2|a it is clear that Optical QSOs (high L AGNs) produce black holes with 
masses Mbh > lO^M© while X-ray AGNs (low L AGNs) produce the lower mass 
ones. Moreover, comparing the BHMF obtained from AGN activity at z = 
with the one cit Zq it can be clearly seen that BH masses are produced during 
the AGN activity: this is the main reason why the results are little sensitive 
to the initial conditions. The model presented in the figure has R2-i=0, i.e. 
it considers only unobscured AGNs. Thus it is not unexpected to derive a low 
pgjj =1.8xlO^M0 Mpc"^. Let's first consider the high L AGNs. It is not yet 
clear if type 2 quasars exist or not or what is their ratio R2-1 relative to type 
1 objects. The comparison between the "observed" local BHMF with the one 
deduced from AGN activity can yield a constraint on R2-1. From Fig. ^3, 
1 < R2-1 < 2 is required to match the local BHMF in absolute level (the slope 
is nicely reproduced). Correspondingly pbr = 3.2-^4.6x IO^Mbh Mpc~^. Let's 
now focus on the low L AGNs. The existence of low L type 2 AGNs is well know 
and the local ratio is R2_i=4 (e.g. Maiolino & Rieke 1995). As shown in Fig. ^, 
with that ratio the local BHMF at low masses (Mbh < IO'^Mq) is overproduced. 
A change in the bolometric correction does not improve the situation because 
the BHMF is just shifted along the M axis. A possible solution to this problem is 
that low luminosity AGNs are not emitting at the Eddington Luminosity (Fig. 
|3|a), a possibility also considered in other works (e.g. Salucci et al. 1999 and 
references therein), and indeed, apart for some "wiggles", the BHMF is now 
well reproduced. One must keep in mind that this might also mean that there 
are some problems with the AGN luminosity functions at low luminosities or 
that the MBH-^Buige relation is different at low masses. 

5. Summary and Conclusions 

We have used the continuity equation to follow the evolution of the BH mass 
function in the assumption that mass accretion onto a BH is the powering mech- 
anism of AGN activity. AGNs accrete at a fraction of their Eddington luminosity 
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and at a starting redshift Zo all BHs were active. By comparing the expected 
BHMF with the one derived from local galaxy luminosity functions we can make 
the following conclusions. (1) The majority of the mass of a local BH was pro- 
duced during AGN activity. (2) Quasar activity was responsible for the growth 
of BHs with M > IO^Mq. (3) The final results are not sensitive to the value of 
the redshift Zo at which all BH were active. (4) The ratio between type 2 and 
type 1 quasars is R2-1 ~ 1 — 2, i.e. ~ 30% — 50% of pbh is made by type 2 
QSOs. (5) The required R2-i._b can be reduced if we assume that L < LEcidington 
or e < 0.1. In any case it is clear that (a) R2-i,_b is much lower than the 4—10 
value required to fit the XRB and (b) no small efficiencies are required to explain 
the local BHMF. (6) R2-i,x > gives a problem at low masses overprcdicting 
the number of relic BHs. This could indicate that either L < LEddington at low 
luminosities or that there might be uncertainties on the BHMF and/or the AGN 
LF at low M and L. (7) No merging of BHs is required by our model. 

The main conclusion of this work is that with "reasonable" assumptions, 
compatible with current knowledge, mass accretion onto massive BH during 
AGN activity can explain the local BH mass function both in shape and nor- 
malization. There is not much room for a large {R » 1) population of high 
luminosity obscured AGNs (i.e. type 2 quasars). 
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